Reducing spurious gravitational radiation in 
binary-black-hole simulations by using conformally 
curved initial data 



Geoffrey Lovelace 

Center for Radiophysics and Space Research, Cornell University, Ithaca, NY, 14853 
Theoretical Astrophysics 130-33, California Institute of Technology, Pasadena, CA 
91125 

E-mail: geof f reySastro . Cornell . edu 

Abstract. At early times in numerical evolutions of binary black holes, current 
simulations contain an initial burst of spurious gravitational radiation (also called 
"junk radiation") which is not astrophysically realistic. The spurious radiation is 
a consequence of how the binary-black-hole initial data are constructed: the initial 
data are typically assumed to be conformally flat. In this paper, I adopt a curved 
conformal metric that is a superposition of two boosted, non-spinning black holes that 
are approximately 15 orbits from merger. I compare junk radiation of the superposed- 
boosted-Schwarzschild (SBS) initial data with the junk of corresponding conformally 
flat, maximally sliced (CFMS) initial data. The SBS junk is smaller in amplitude than 
the CFMS junk, with the junk's leading-order spectral modes typically being reduced 
by a factor of order two or more. 



PACS numbers: 04.25.D-,04.25.dg,04.20.Ex,02.70.Hm 
1. Introduction 

One of the most important sources of gravitational waves for LIGO [1] is the inspiral 
and merger of two black holes. LIGO has reached its design sensitivity and can detect 
stellar-mass binary-black-hole mergers as distant as about 100 megaparsecs [2]. 

The gravitational waveform of a binary-black-hole merger cannot be computed 
using analytic techniques but must be obtained by solving the Einstein equations 
numerically. Currently, simulations of binary black holes are based on splitting the 
four-dimensional spacetime into a series of three-dimensional spatial slices; to start an 
evolution, one must construct initial data for the first slice. This initial data must i) 
represent the desired physical situation (i.e., two black holes about to merge), and ii) 
satisfy the vacuum Einstein constraint equations: 

Gnj = 0, (1) 
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where the subscript n refers to components normal to the initial slice and the subscript 
j refers to components tangent to the slice. The normal-normal and normal-spatial 
equations are called the Hamiltonian and momentum constraints, respectively. The 
vacuum Einstein evolution equations, Gij = 0, are solved to step from the initial slice 
to subsequent slices. 

There are several methods that generate constraint-satisfying initial data (for a 
review, see, e.g., |3]), including Bowen-York puncture data [HIS] and quasiequilibrium 
excision data [SI [71 El El [IQ] . However, these methods generally assume that the initial 
spatial metric Qij is conformally flat: 

= i^'h. (2) 

where fij is the metric of fiat space. This simplifying assumption causes spurious 
gravitational radiation (also called "junk radiation") to be present in the early phases 
of the simulation. Specifically, it is known that a stationary, isolated black hole with 
linear [11] or angular [12] momentum cannot be sliced so that the spatial metric is 
conformally fiat. Attempting to construct constraint-satisfying, conformally fiat initial 
data for boosted or spinning black holes yields holes that are not in equilibrium but 
are unphysically perturbed. As they relax to an equilibrium configuration, gravitational 
waves are emitted. 

In general, the black holes in a binary have both linear and angular momentum; 
therefore, binary-black-hole simulations using conformally fiat initial data will also 
contain junk radiation. Before one can extract the physically relevant gravitational 
wave signal, one must first evolve the unphysical system until the spurious waves have 
left the computational domain; thus junk radiation adds to the already considerable 
computational expense of the simulation. But besides the additional cost, junk 
radiation causes several other undesirable effects. Spurious gravitational radiation can 
unrealistically shorten the time until the black holes merge [13]. Junk radiation also 
makes accurate comparison with post-Newtonian waveforms more difficult: it reduces 
the accuracy of the simulations and postpones the starting time at which the post- 
Newtonian comparison can begin [TH [15] . For simulations of binary black holes which 
recoil, the spurious gravitational radiation and astrophysically-realistic gravitational 
waves both carry linear momentum in the direction of the kick; consequently, before 
integrating the radiated momentum fiux to obtain the kick velocity, one must first wait 
for the junk radiation to leave the computational grid [HI [T7] . 

To minimize these effects, it is desirable to reduce the spurious radiation as much as 
possible. Because conformal flatness is known to contribute to the junk radiation, one 
common approach is to choose a curved conformal metric; suitable choices lead to smaller 
amounts of junk radiation. (Note, however, that even conformally flat initial data can be 
constructed such that less spurious radiation is emitted than for standard Bowen-York 
puncture data [IH]-) In Ref. [I9], the authors (building on the work in Refs. [201 [211 [22]) 
have constructed and evolved conformally curved initial data to reduce the amount of 
junk radiation in head-on mergers of spinning holes; by using a conformal metric that is 
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a superposition of two spinning black holes, they find that the leading order, quadrupole 
spherical-harmonic mode of the spurious radiation is reduced. Furthermore, Ref. [23] 
proposes using conformal metrics made by asymptotically matching post-Newtonian and 
perturbed-black-hole metrics, and initial data based on post-Newtonian free data have 
been constructed [211 [25], [26], [27|; however, to the best of my knowledge, no numerical 
simulations evolving these initial data have been published to date. 

In this paper, I construct conformally curved initial data for two nonspinning black 
holes of equal mass in orbits with low eccentricity. In particular, I use a conformal 
metric [Eq. (1271) ] that is a superposition of two boosted Schwarzschild black holes. To 
construct binary-black-hole initial data with nearly-extremal spins, Ref. [28] adopts a 
conformal metric similar to the one used in this paper but superposes two boosted, 
spinning, Kerr-Schild black holes instead of two nonspinning, boosted, Schwarzschild 
black holes. Note that the superposed metrics in this paper and in Ref. [28| are similar 
to the superposed-Kerr-Schild conformal metrics used in Refs. [291 [30l [31] ; however, the 
conformal metrics used here and in Ref. [28] are flat everywhere except near each black 
hole. 

After choosing superposed-boosted-Schwarzschild free data, I then combine them 
with quasiequilibrium boundary conditions developed by Cook [8], Cook and Pfeiffer [9], 
and Caudill et. al. [10]. I solve the constraint equations using the Caltech- Cornell 
pseudospectral elliptic solver [32]. After reducing the eccentricity of the holes' orbits 
by using the technique of Pfeiffer et. al. [33] (which was extended to the conformally 
curved case in Ref. [2H]), I evolve the holes using the Caltech-Cornell code [3ll 135] . 

The remainder of this paper is organized as follows. In Sec. 121 I summarize the 
formalism that I use to solve the initial value problem. In Sec. IHl I describe how 
to construct initial data whose conformal metric is a superposition of two boosted 
Schwarzschild black holes. In Sec. 14.11 I choose a conformally flat, maximally sliced 
(CFMS) initial data set and a superposed-boosted-Schwarzschild (SBS) initial data set 
that is physically comparable. In particular, I use the eccentricity-reduction technique 
of [33] (originally developed under the assumption of conformal flatness) so that sets 
CFMS and SBS both have very little orbital eccentricity. Evolutions of the conformally 
curved data are discussed in Sec. 14.21 and the junk radiation of sets CFMS and SBS are 
compared in Sec. 14.31 A brief conclusion is made in Sec. IS 

2. The initial value problem 

2.1. The constraint equations 

To construct constraint-satisfying initial data, I begin with the usual 3+1 split, in which 
the four- dimensional spacetime, with metric g^y, is split into a series of three-dimensional 
spatial slices with spatial metric gij. The spacetime metric g^y is related to the spatial 
metric gij., the lapse a, and the shift by 

ds^ = g^ydxi'dx" = -a^dt^ + gij{dx' + (3'dt) (dx^ + f3^dt). (3) 
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Here and throughout the rest of this paper, the Einstein summation convention is 
assumed. Greek indices refer to spacetime coordinates and are raised and lowered with 
the spacetime metric g^v and its inverse. Latin indices refer to spatial coordinates of a 
t = const slice and are raised and lowered with the spatial metric gij and its inverse. 

On the initial [t = 0) slice, the initial data must specify gtj and the extrinsic 
curvature Kij, which is essentially the rate of change of gij in the normal direction. The 
extrinsic curvature is related to the time derivative of the metric dtgij and to the lapse 
and shift by 

dtg,, = -2aKij + 2V^i(3,y (4) 

The initial values of gij and Kij must be chosen so that i) the solution contains the 
desired physical content, and ii) the constraint equations ([H) are satisfied. 

A systematic way to solve these equations is given by the extended conformal thin 
sandwich (XCTS) formalism |36| [37] . In this formalism, one expands gij and Kij as 
follows: 



9ij = gij, 

Ki, = A, + ^-g,,K. (5) 

Then, one chooses the conformal metric cjij, the trace of the extrinsic curvature K, and 
the time derivatives of both, Uij := dtgij and dfK. The constraint equations ([1]) are then 
reduced to elliptic equations for the conformal factor ip and the shift A fifth elliptic 
equation for determines the lapse; it is not a constraint, but appears because the free 
data include dtK instead of a := i^i'^a. (Alternatively, one could use the "standard" 
conformal thin sandwich equations [36l [37], in which the free data are gij, Uij, K, and 
a.) 

Together, these equations form a second-order, nonlinear, coupled elliptic system 
called the extended conformal thin sandwich (XCTS) equations, which are (e.g., Eq. (8) 
of [33]): 

- - l^V' + li^-'A^'Aj = 0, (6a) 



V^aij) - (aij) \^ + ^K^ij^ + l^"^A'^AJ = -^j^dtK - p'd.K). (6c) 
.olZ o 

Here V is the gradient with respect to gij, the "longitudinal operator" L is twice the 
symmetric, trace-free gradient (i.e., the "shear") with respect to gij, i.e., 

(L^),, := V,V, + VjV, - h^.VkV", (7) 



and 



A'^ = ij^^A'^ = -— - {hpy^ - u'^ . (8) 
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The initial value problem now amounts to i) choosing the free data ( cjij, Uij, K, and 
dfK), ii) choosing boundary conditions for ip, aip, and P\ and iii) solving Eqs. fl6aj) - fl6cj) 
for ijj, aip, and /3\ Most of these have preferred choices, motivated by the requirement 
that, in the comoving coordinates, the initial data contain two black holes at rest. These 
quasiequilibrium conditions will be discussed in the next subsection. The remaining 
quantities will be dealt with in Sec. [31 

2.2. Quasiequilibrium free data and boundary conditions 

In the XCTS formalism described in the previous subsection, the physical content 
of the data is selected by the choice of both the free data {cjij, K, and their time 
derivatives) and by the boundary conditions. We wish to make choices that represent 
the physical situation of two (otherwise isolated) black holes orbiting each other. In 
the quasiequilibrium method [6], O [H O [TOl [33] used in this paper, there are preferred 
choices for many of the free data and boundary conditions. 

2.2.1. Free data In quasiequilibrium initial data, the coordinates are required to 
(initially) be comoving with the black holes. If the holes are also in equilibrium, time 
derivatives in the comoving frame should initially be small. Quasiequilibrium initial 
data therefore choose 



The remaining free data, gij and K, can be chosen freely. In Sec. [3l I make particular 
choices for gij and K. 

2.2.2. Outer boundary conditions The computational domain can be represented by 
only a finite number of gridpoints, so it necessarily will have an outer boundary B, which 
here is taken to be a coordinate sphere whose radius R is so much larger than all other 
length scales that it is effectively "infinitely far away." (In practice, the outer boundary 
is roughly 10^ times larger than the size of each black hole.) The physical requirement 
that the binary is isolated (i.e., that the spacetime is asymptotically flat) corresponds 
to the conditions 



provided that gij is asymptotically flat. 

The outer boundary condition on the shift is set by the requirement that the 
coordinates are initially comoving with the black holes. Therefore, in the asymptotically 
flat region — and in particular, on B — the coordinates will not be inertial; instead, they 
will rotate (due to the orbital motion) and contract (due to the holes' inspiral). That 
is, if r is a coordinate radius measured from the system's center of energy, and if r* is a 



Uij = 0, 

m = 0. 



(9) 
(10) 



ifj = 1 on B, 
a = 1 on i5. 



(11) 
(12) 
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radial position vector in the asymptotically fiat region, then 

(5' = (flo X r)^ + —r' on 8. (13) 

ro 

Here ro = do/2, where do is the initial coordinate separation of the holes. The precise 
values of Qq and Vr will be set so that the holes' subsequent trajectories are not eccentric 
(Sec. ED. 

2.2.3. Inner boundary conditions The singularities of each black hole are excised from 
the computational domain. The excision surface S is chosen to be the apparent horizons 
Ti of the two holes (labeled "A" and "B"), i.e., S = HaIJ'Hb- This requirement leads 
to a boundary condition on the conformal factor (Eqs. (28) and (48) of [9]): 



S''^k^/J = -^s's' {hf3)ij~u,j -I^h'^WiSj + ^KiP' on S (14) 



6a 



4 ' ^ 6 



where s* is an outward-point ingji unit normal vector on S, and hij = gij — SiSj = 
{9ij ~ SiSj) = ip'^hij is the induced metric on S. 
When the initial data are evolved in the comoving system, the apparent horizon 
(itself in equilibrium) should remain at rest. This leads to the following boundary 
condition on the shift (Eqs. (36) and (50) of [9]): 

13' = as' + on S. (15) 

Here is a parameter that determines the amount of spin on the hole in addition 
to corotation, and C,' is a conformal Killing vector within S. (If the holes are to have 
different spins, different values of Qr are used on TCa and Hb-) The (/c)**^ component of 
the quasilocal spin [381 EHl HQ] of each hole can be written as (e.g., Eq. (37) of [TO] ) 

(fc) 

o" Jn 

where Ti is either Ha otT-Cb, and ^^^-j is a Killing vector on Ti. that defines rotations about 
the k axis. In general, there are not necessarily any Killing vectors on 7i; in such cases, 
the blacli hole spin can still be measured using Eq. (fT6|) but with taken to be an 
approximate Killing vector [HI |12] . In this paper, I consider only non- rotating binaries, 



^ £ (K,, - g,,K) q,/S\ (16) 



in which Qj. is selected so that the approximate-Killing-vector spin is very close to 
zero. (The specific method that I use to measure the black-hole spins is described in 
Appendix A of Ref. [28].) 

The inner condition on a is a gauge choice [9]; i.e., it does not affect the physical 
content of the initial data. The particular choices used in this paper are discussed in 
Sec. [31 

I Here "outward-pointing" points away from the black hole, toward infinity. 
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3. Nonspinning, non-eccentric binary-black-hole initial data 

3.1. Conformally flat data 

Initial data for binary-black-holes are typically assumed to be conformally flat: 

9ij = fir (17) 

In Sec. 14.11 the conformally flat, maximally sliced data set CFMS (which is "30c" in 
Ref. [H]) chooses 

K = 0, (18) 

dr{a^) = on 5. (19) 

3.2. Conformally curved data 

To reduce the amount of junk radiation, I adopt a curved conformal metric. In this 
section, I build up a suitable conformal metric for two non-spinning black holes in an 
initially circular orbit. 

3.2.1. Schwarzschild with maximal slicing The Schwarzschild metric can be split into 
maximal slices (i.e., slices with K = 0.) With maximal slicing, the Schwarzschild 
spacetime is (e.g., Eq. (52) of [9]): 

g^jdx'dx^ = a^^dR^ + R^dn\ 



= (20) 

where R is the Schwarzschild areal radial coordinate and e^j is a unit vector normal to 
constant-i? surfaces. The choice of C specifies which maximal slicing is used (and thence 
the coordinate radius of the horizon). To facilitate comparison with [33], in this paper 
I choose C = 1.737, which implies the horizon radius is at rg^c = 0.8595, and Ms = 1. 

3.2.2. Metric of a boosted Schwarzschild hole To construct the metric of a boosted, 
Schwarzschild hole, I begin with the Schwarzschild metric [Eq. (I20l)]. Then, I apply 
the following two coordinate transformations: 

(i) First, I apply a radial coordinate transformation to make the spatial metric 
conformally flat. Any spherically-symmetric metric can be made conformally flat 
by i) writing it as 

g.jdx'dx^ = f{R)dR^ + R^d^'^, (21) 

where R is an areal radial coordinate, and then ii) making a radial coordinate 
transformation r = r{R) so that the metric is conformally flat: 

f{R)dR^ + R^dn'' = {dr^ + r^dn^) (22) 
^ = i?/r and |^ = ^V^- (23) 
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For a given f{R), determining r{R) amounts to solving a first-order ordinary 
differential equation. For the choice f{R) = [cf. Eq. f l20|) ]. the analytic solution 
to the ODE is unknown, so the ODE ( l23l) is solved numerically. 

(ii) Once the spatial metric is conformally flat, I give the Schwarzschild hole a velocity 
V in the y direction by performing a Lorentz boost. The resulting spacetime metric 
is g^r\v). 

3.2.3. Superposing two boosted, non-spinning holes Suppose two nonspinning black 
holes, "A" and "B," are initially separated by coordinate distance do along the x axis. 
Suppose they are in a circular orbit about the z axis with speed v = QQdo/2. 

I seek free data that accurately describes the boosted black holes. The simplest 
choice is to merely superpose two boosted-Schwarzschild black-holes: 

~9^, = k + (4 - k) + (4 - ' (24) 
K =K'^ + K^. (25) 

For hole A, gf^ and are obtained by translating g^°"^^{v) so that the hole is 
centered about (x, y, z) = {do/2, 0, 0). Likewise, for hole B, gfj and are obtained by 
translating so that the hole is centered about {x,y,z) = {—do/2,0,0). Note 

that as V ^ 0, cjij — ^ "^ofij K ^ 0; after absorbing ipo into the conformal factor ip, 
the V = superposed-boosted-Schwarzschild data in fact reduces to the well-examined 
(conformally flat) quasiequilibrium data by Cook and Pfeiffer [9]. 

I used the Caltech-Cornell pseudospectral code [32] to solve the XCTS equations 
with the free data and boundary conditions described previously. When the conformal 
metric gij is given by Eq. (!24|) . I find that the elliptic solver does not converge [Fig. [1]. 

The source of the difficulty can be seen by inserting the outer boundary conditions 
ffTTj) . f|T2|) . and f|T3l) into the Hamiltonian constraint [the first of Eq. ([1])], which gives 
an equation of the form 

- V'V' - [L (rJo X r)]'^' [L {Qo x r)],, + ■ ■ ■ = 0. (26) 

The source term shown here vanishes when the metric is conformally fiat, because 
(fio X r)* is a conformal Killing vector of fiat space. But when the conformal metric is 
the naive superposition (l2lll . the term [L (fio x r)]^'^ contains a spherically-symmetric 
part that decays only as 0{l/r) [Fig. [2]. It follows that const x logr as r ^ oo, 
but this is incompatible with the requirement that ■?/; ^ 1 as r — * oo. 

3.2.4. Scaling the non-flat terms by Gaussians The inconsistency described in the 
previous subsection can be avoided by requiring that the non-fiat terms fall off 
sufficiently quickly far from the holes. This can be accomplished by scaling the non-fiat 
terms by Gaussians: 

~m = ■■= k + e-''^^-' {gf, - k) + e-y-' (4 - /^,) , (27) 
K = i^s^s _ e-r\/u,^K'^ + e-"s/-'j^5. (28) 
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Constraints vs resolution (naive superposition) 




Figure 1. Color online. When the conformal metric cjij is given by Eq. p4|) . the 
constraints do not decrease exponentially with increasing resolution (indicated here by 
the number of gridpoints N on the computational grid). 

Here "SBS" stands for "superposed, boosted Schwarzschild," and and are the 
coordinate distances from the centers of holes A and B, respectively. The width w can 
be adjusted to minimize the junk radiation; further investigation is needed to find an 
optimal choice for w. The Gaussians used here may seem similar to the attenuation 
functions of [31]; however, the attenuation functions here approach zero at large radii, 
whereas those of [31] approach unity. In Sec. 14.11 data set SBS uses w = 20rexc- With 
this choice, the elliptic solver converges, and the constraints decay exponentially with 
resolution, as expected [Fig. [3]. 

Finally, as noted in Sec. 12.2.31 the inner boundary condition on the lapse a is a 
gauge choice; in set SBS, the condition is 

= 1 + e""^/"' (a^ - l) + 6"^^/"'' (a^ - l) on S. (29) 

Here is obtained by translating g^°°^^{v) so that hole A is centered about (x, y, z) = 
{do/2, 0, 0); likewise, is obtained by translating g^^°^^{—v) so that hole B is centered 
about {x,y,z) = {—do/2,0,0). 
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Spherically symmetric part of shear of {Q. x r) vs radius 
Confomal metric = Naive superposition 




Coordinate radius 



Figure 2. Color online. Magnitude of selected components of [L (ilo x r)]*-' ; the 
plot shows the absolute value of the angular average as a function of radius for each 
Cartesian component. The component [L {Q,q x r)]^** includes a spherically-symmetric 
term that decays only as 1 /r when the conformal metric is given by Eq. (j24p ; this term 
causes the conformal factor to diverge logarithmically as r ^ oo. The other nonzero 
components decay as l/r^ or faster. 



Quantity 


CFMS 


SBS 




0.016708 


0.016405 


Vr 


-4.3 X 10-^ 


-3.5 X 10"^ 


Mi,, 


1.039 


1.027 


M := 2Mi„ 


2.078 


2.044 


Madm 


2.062 


2.039 


la^/Mi^J and ja^/M^J 


< 2 X 10-^5 


< 2 X 10"^° 




< 6 X 10-s 


< 5 X 10^6 


do/M 


14.44 


14.60 


So/M 


17.37 


17.51 



Table 1. A comparison of the two initial data sets presented in this paper. Set CFMS 
is the conformally flat initial data set "30c" in Ref . [2] , and set SBS uses a conformal 
metric that is a superposition of two boosted Schwarzschild black holes. The initial 
data sets describe physically comparable situations: the masses and separations agree 
to within about 1%, and the frequencies agree to within about 2%. The radial velocities 
are comparable [and are chosen so that the eccentricity is small (Fig. [4])] , and the spins 
of the holes are close to zero in both cases. 
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Constraints vs resolution (SBS) 




Figure 3. Color online. Constraints as a function of resolution (indicated by the 
number of gridpoints N) when the conformal metric is gf^^^ which is curved near the 
black holes but fiat far away. The constraints decrease exponentially with increasing 
resolution, as expected. 



4. Comparing the junk radiation of conformally flat and 
superposed-boosted-Schwarzschild initial data 

4.1. Initial data sets 

In this section, I compare the amount of junk radiation during evolutions of two initial 
data sets, both of which correspond to equal-mass, non-spinning binary black holes in 
nearly-circular obits; the holes are about 15 orbits away from merger. The data sets are 
i) a conformally flat, maximally sliced set (CFMS) which is identical to the initial data 
set "30c" that is constructed and evolved in Ref [H], and ii) a superposed-boosted- 
Schwarzschild set (SBS) which is constructed using the method described in Sec. |3l 

Table [H compares some physical properties of the two data sets. The table lists 
three measures of mass: i) the irreducible mass Mi„ = ^J~Ap^[JiQtx , where Aah is the 
area of one hole's apparent horizon, ii) the sum M of the holes' irreducible masses, 
and iii) the Arnowitt-Deser-Misner (ADM) Mass, a measure of the total energy in the 
system, defined by the following surface integral at spatial infinity: 

MADM = Y^y" {d,g,,-dig,,)d''S\ (30) 

In practice, the integral for Madm is evaluated on the outer boundary of the initial-data 
grid (i.e., on i3, which is a sphere whose radius is approximately lO^Mjrr). As shown 
in Table [H both Mirr and Madm agree within about 1%. Note that in the rest of this 
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-- SBS: eccentricity reduction iter. 1 
— SBS: eccentricity reduction iter. 2 



200 400 600 800 1000 1200 1400 

Time / M 

Figure 4. Color online. Reducing the eccentricity in the SBS initial data set. The 
proper separation s vs. time is shown for the conformally flat, maximally sliced 
(CFMS) data set (dashed red line). The corresponding SBS data set with the same 
choices for fig and Vr is much more eccentric (dotted black line), but two iterations 
of the algorithm described in Sec. 4 of Ref. [33] greatly reduce the orbital eccentricity 
(dashed and solid black lines). All evolutions shown in this plot were performed at 
resolutions of 46.8'^ gridpoints. 



section, following Ref. [H], I will typically express quantities with dimension in terms 
of the mass M. 

The spins of the holes are computed using the quasilocal approximate-Killing- vector 
spin [Eq. f fT6l) and the surrounding discussion]. Tabled] list the quasilocal spin for both 
the CFMS and SBS data sets; in both cases, the spin is close to zero. 

The initial coordinate separation of the holes d^ is given for sets CFMS and SBS 
in Table [U Also given is the initial "proper separation" So, which is defined in terms of 
the following line integral along the x-axis of the comoving frame: 

/l-d/2—rexc 
ds= dx^/g^ (31) 

J ~d/2+rexc 

Here the limits of integration are the coordinate locations where the holes' apparent 
horizons intersect the x axis, and the radius Ve^c is the coordinate radius of the apparent 
horizon. Note that the "proper separation" is coordinate-dependent (since the integral 
is taken along the x-axis, rather than along a geodesic of extremal length) and slicing- 
dependent (since the distance is measured within the spatial slice). 

The initial orbital frequency Qq and radial velocity Vr listed in Table [1] are chosen 
to reduce the holes' orbital eccentricity. The eccentricity reduction for set CFMS is 
described in Sec. IIA of Ref. [H]. To reduce the eccentricity of the SBS data, 1 initially 
guess that Qq and Vr have the same coordinate values as the CFMS values; then, 1 
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Figure 5. Color online. The constraint violation plotted as a function of time for 
evolutions of conformally flat, maximally sliced (CFMS, left panel) and superposed- 
boosted-Schwarzschild (SBS, right panel) initial data. The constraints are shown for 
three different spatial resolutions (labeled by their corresponding number of gridpoints 
N). 



tune them using the iterative scheme described in Sec. 4 of [33]. Each iteration would 
completely remove the eccentricity if the binary were Newtonian; in the relativistic case, 
successive iterations converge to non-eccentric orbits. Figure H] illustrates the eccentricity 
reduction by showing the proper separation as a function of time for set CFMS and also 
for several iterations of SBS data. The set "SBS: eccentricity reduction iter. 2" is the 
set SBS described in Tabled] and used in the remainder of this paper. 

4-2. Evolutions 

The CFMS and SBS data sets were evolved using the Caltech-Cornell pseudospectral 
evolution code SpEC [35]. 1 used the same evolution methods, equations, and boundary 
conditions as those described in Ref . [11] . In particular, 1 use the first-order generalized- 
harmonic system with constraint damping that is derived in Ref. [31]. The outer 
boundary conditions, derived in Ref. [3l] and augmented with improved conditions on 
the gauge fields in Refs. [131 HI] , preserve the constraints and enforce a "no- incoming- 
radiation" requirement by freezing the Newman-Penrose scalar ^>q. 

The evolution grid's outer boundary is at approximately QOOMirr. The excision 
boundaries are slightly inside the apparent horizons; this is accomplished by 
extrapolating the initial data from S to points slightly inside S. (For set CFMS, 
the evolution-grid excision spheres are at radius rg^, = Q.QTrexc] to accommodate 
nonspherical horizons, the evolution-grid excision spheres are at rev = 0.93rexc for set 
SBS.) Aside from differences in extrapolation, the evolutions of sets CFMS and SBS 
used identical computational domains and identical numerical resolutions. 

Each evolution of sets CFMS and SBS were performed at three different numerical 
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Figure 6. Color online. A comparison of the dominant (quadrupolar) modes of 
tlie junlt radiation for conformally flat (CFMS, left panel) and superposed-boosted- 
Schwarzschild (SBS, right panel) initial data. The (£,m) = (2,0) mode (red) and 
{£,m) = (2,2) mode (black) of |^'4| are shown at at three different spatial resolutions: 
46. 8^ (dotted Uncs), 51. 9^ (dashed hncs), and 57. 0^ (sohd lines). The waves are 
extracted on a coordinate sphere of radius 265A/. The blue, dashed horizontal line 
is the maximum amplitude of the SBS junk. 



resolutions, corresponding to 46.8^, 51.9'^, and 57.0'^ gridpoints, respectively. Figure [5] 
compares the constraint violations (which are computed in the same way as in footnote 
8 of Ref. [28]) at each resolution. At late times (> 500M), after the initial burst of junk 
radiation has left the computational domain, the constraints converge exponentially, 
as expected — but while the junk radiation is on the grid, the constraints do not 
appear to decrease significantly between the = 51.9'^ and = 57.0^ resolutions. 
This poor convergence is not surprising; because junk radiation typically has a much 
higher frequency than the physical gravitational waves, the spurious radiation requires 
much higher numerical resolutions to be accurately resolved. In practice, such high 
resolutions are never used, because the added computational cost (which would likely 
be prohibitively expensive) is not necessary to adequately resolve the astrophysically- 
realistic content of the simulation. 

Because my motivation is to reduce the observed amount of spurious radiation, 
not to study its physical content, in this paper I do not attempt to better resolve the 
junk radiation. Instead, in the next subsection, I simply compare the amounts of junk 
radiation that can be seen for the three resolutions shown in Fig. [51 

4-3. Junk radiation comparison 

The gravitational waves are extracted from the CFMS and SBS evolutions at the same 
coordinate radius rextr- Specifically, the simulation computes the Newman- Penrose 
scalar \l/4, which at large radii is related to the outgoing +- polarized and x -polarized 
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Figure 7. Color online. A comparison of tlie gravitational waves extracted at 
coordinate radius roxtr = 265M for the conformally flat, maximally sliced (CFMS, 
left panel) and superposed-boosted-Schwarzschild (SBS, right panel) initial data. The 
(. = 2 spherical-harmonic modes are shown in black, and higher modes are shown in 
red. All of the modes of |^'4| shown here were computed with a spatial resolution of 
N ~ 57.0"^ gridpoints. The blue, dashed, horizontal line is the maximum amplitude of 
the SBS junk. 



^4 = 3T^/^+-^;^/^x. (32) 



gravitational- wave amplitudes by 

d^ rf2 

The scalar \E'4 is evaluated on a sphere of radius rgxtr = 265M and then expanded in spin- 
weighted-spherical-harmonic modes (For further details on the wave-extraction 
method used here, see Sec. 5.3 of Ref. [35].) 

At early times, the waveform consists of spurious gravitational waves; they are 
recognizable as such by their frequencies, which are much higher than the dominant 
frequencies / ~ f2o/7r of the physical, quadrupolar gravitational waves. Figure [6] plots 
dominant, quadrupolar components of Mr extr l^^h they have a frequency / ~ (15M)~^, 
which is significantly higher than the physical frequency / ~ VLq/ti ~ (200M)~^. 
Figure [7] compares all of the spherical harmonic modes of Mr extr |^^| from £ = 2 through 
1 = 7 with m > 0. 

Following the initial high-frequency spurious radiation, the astrophysical waves are 
dominated by the (£, m) = (2, 2) mode, with Mr extr 1^22! ^ dimensionless slowly-growing 
amplitude that is essentially constant during the time shown in Fig. [71 during this 
interval, this amplitude is significantly smaller than the peaks of the junk radiation. 

As a simple measure of the amount of junk radiation present in each {i,m) mode. 
Fig. [8] shows the maximum value of Mr extr l^^ml modes from i = 2 through 

i = 5. Only modes with m > are shown. The numerical values of the maxima 
vary with increasing resolution; this is another indication that the junk radiation is 
underresolved. However, for each resolution, the SBS maximum is smaller than the 
corresponding CFMS maximum. For most of the modes with large contributions to the 
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Figure 8. Color online. A comparison of the maximum amplitude of tlic spherical 
harmonic modes of |5'^,„| during the time interval shown in Figs. [BHT] The [i, m = 0) 
modes are labeled by £, and the m > 0) modes are indicated by small tick 
marks (with m increasing from left to right). The three panels correspond to three 
different spatial resolutions, with resolution increasing from left to right. The waves 
are extracted at a radius of rcxtr — 265M. At each resolution, the nonvanishing 
conformally flat, maximally sliced (CFMS, red, thin line) modes are larger than the 
corresponding superposed-boosted-Schwarzschild (SBS, thick, black line) mode. Only 
the modes up through I = b with m > are shown. 

junk radiation — including the {i,m) = (2,0) mode, which is the leading contributor 
to the junk — the improvement is about a factor of two or more; however, for the 
{i,m) = (2,2) mode, the improvement is much more modest: the SBS maximum is 
about 75% of the corresponding CFMS maximum. The smaller improvement might be 
because the SBS initial data makes no attempt to model the initial gravitational-wave 
content; if so, then the spurious gravitational radiation might be able to be reduced 
even more by adding post-Newtonian terms to gij and K (such as, e.g., terms based on 
the free data used in Ref. [27]). 

5. Conclusion 

The junk radiation in binary-black-hole simulations can be significantly reduced by using 
superposed-boosted-Schwarzschild initial data instead of conformally flat data. For the 
case of two non-spinning black holes initially 15 orbits from merger, the amplitude of 
the junk gravitational waves decrease, with most spherical-harmonic modes, including 
the leading contributors to the junk radiation, decreasing by a factor of order 2 or 
more. However, a significant amount of junk radiation is still present in evolutions of 
the superposed-boosted-Schwarzschild data, which also does not attempt to account for 
the initial gravitational- wave content. Besides attempting to incorporate gravitational 
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waves into the initial data, future studies could investigate the amount of spurious 
radiation in the superposed-Kerr-Schild data in Ref. [28], since one might speculate that 
lower amounts of spurious radiation play a role in letting the superposed-Kerr-Schild 
data lead to evolutions with larger spins than are possible with Bowen-York initial data. 
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